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Abstract 

Static and spherically symmetric perfect fluid solutions of Einstein's 
field equations with cosmological constant are analysed. After showing 
existence and uniqueness of a regular solution at the centre the extension 
of this solution is discussed. Then the existence of global solutions with 
given equation of state and cosmological constant bounded by 4npb, where 
pb is the boundary density (given by the equation of state) of the perfect 
fluid ball, is proved. 
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1 Introduction 



This paper analyses static, spherically symmetric perfect fluid solutions to Ein- 
stein's field equations with cosmological constant for a given monotonic equation 
of state p = p(P). The choice of central pressure (central energy density) and 
cosmological constant uniquely determines the pressure function. It is the aim 
of this work to extend the results presented in ^2] to include the cosmological 
constant. 

Existence and uniqueness of a global solution with given equation of state 
can be proved for cosmological constants satisfying A < Anpb, which is given by 
the equation of state since pb = p(P = 0), following the line of argument of j!2| . 

The existence of such global solutions is quite important because assuming 
their existence, Buchdahl ;5 showed that the total mass of a fluid ball is bounded 
by its radius. He showed the strict inequality M < (4/9)i?, which holds for fluid 
balls in which the density does not increase outwards. It implies that radii of 
fluid balls are always larger than the black-hole event horizon. It is expected 
that perfect fluid balls with cosmological constant are larger than the black hole 
event horizon but still be smaller than the cosmological event horizon. 

This paper is organised as follows: In section[5]Einstein's field equations and 
the Buchdahl variables are presented, in the beginning of section [21 uniqueness 
and existence of a regular solution at the centre is shown. In theorem El the 
extension of the solution is discussed. For fluid balls the Buchdahl inequality 
with cosmological term is obtained in theorem^ Solutions without singularities 
are constructed in section [5J Finally remarks on the finitcness of solutions are 
given and moreover some conclusions are presented. 

2 Field equations and Buchdahl variables 

The most general static and spherically symmetric 

ds 2 = -e^dt 2 + e a{r) dr 2 + r 2 {d6 2 + sin 2 9d<j) 2 ) , (1) 

in Einstein's theory of gravity yields three independent field equations with 
cosmological constant 

2- e ^)A. ( r _ re ~<r)\ „ Ae *(r) = 87rp ( r ) e "(r) ( 2 ) 
r z dr \ I 

-4 (l + rv'{f) - e a(r) ) + Ae a M = 87rP(r)e a M , (3) 
-^P-(P(r)+p(r))=P'(r), (4) 

for an isotropic perfect fluid. Note that one can either use the three field equa- 
tions which imply conservation of energy-momentum or two of the field equa- 
tions together with the conservation equation which we do. 

One function can be chosen freely since these are three independent ordinary 
differential equations for four unknown functions. The most physical assumption 
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is to prescribe an equation of state p = p(P). Integration of {2J gives 

e~ a[r) = 1 - 2w(r)r 2 - -r 2 , (5) 
3 

where the constant of integration is pnt to zero demanding regularity at the 
centre and where w(r) is the mean density up to r defined by 

w ( r ) = — V" ' TO ( r ) = 4tt / s 2 p(s)ds . (6) 



Eliminating the function v'(r) from the field equations yields the Tolman- 
Oppenheimer-Volkoff HUHniini equation (TOV-A) 



^ (47rP(r ) + w(r) ~ |) (P( r) + p{r)) 
1 - 2w(r)r 2 - | 

With given equation of state the conservation equation leads to 



where P c is the central pressure. With m(r), equations 10 and JJJ arc forming 
an integro-differential system for p(r) and P(r). However, differentiating the 
mean density w{r) with respect to r implies 

w '( r ) = - Unp(P(r)) - 3w(r)) . (9) 
r 

Therefore for a certain p = p(P) , equations Q and © are forming a system of 
first order differential equations in P(r) and w(r). 

To extract the TOV-A equation from Einstein's field equations the metric 
function v{r) was eliminated. On the other hand, one may eliminate the pressure 
for which Buchdahl jS] introduced new variables 

y 2 = e-°W = i _ 2w(r)r 2 - ^r 2 , (10) 

C = e v/2 , (11) 
x = r 



(12) 



which we have supplemented by a cosmological constant. The second field 
equation @ written in those variables leads to 

2 C 

HttP- -A = -2to. (13) 

After differentiating this with respect to x and using the conservation equation 
fifl the pressure can be eliminated. After some algebra one arrives at 

(vU, x -l^ = o, (M) 
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which is surprisingly similar to the square of the Weyl tensor 

n n abcd 64 V 2 f t > \ ,lw,x(\ 2 2 f-i r\ 

CabcdC = \{y(,x),x + — ) X ■ I 15 ) 

This implies that the constant density solutions are conformally flat since 
w, x = and therefore with i|14[l gives (yC.x).x = 0. 



3 Existence of unique regular solutions 

Since the field equation for a static and spherically symmetric perfect fluid with 
given equation of state reduce to a system of singular first order differential 
equations it is our first aim to apply theorem 1 of |12j to the resulting system 
with cosmological constant. 

Theorem 1 (Rendall and Schmidt, 1991) Let V be a finite dimensional 
real vector space, N : V — > V a linear mapping, G : V x I —> V a C°° mapping 
and g : I — > V a smooth mapping, where I is an open interval in R containing 
zero. Consider the equation 

sf s +Nf = sG(s,f(s))=g(s), (16) 

for a function f defined on a neighbourhood of and I and taking values in V. 
Suppose that each eigenvalue of N has a positive real part. Then there exists an 
open Interval J with € J C / and a unique bounded C 1 function f on J \ 
satisfying \ lh}) . Moreover f extends to a C°° solution of on J if N,G and 
g depend smoothly on a parameter z and the eigenvalues of N are distinct then 
the solution also depends smoothly on z. 

Note that 10 is singular at the centre whereas (7J is not, however using 
p = p c + xpi , where p c is the central density, we find for a given equation of 
state 



.r 



dw 3 
dx + 2 



w = 2irp c + 2-Kxpi , (17) 



dpi , (dP\ 1 l(4nP + w-±){P + p c + x Pl ) 

x ~r- + = - \-r) o : — ir a • ( 18 ) 



dx \ dp J 2 1 - 2wx - jx 

With the help of the corresponding pressure relation P = P c + xP\(pi) and by 
noting that 

A V 1 ( A 

1 — 2wx — — xj = 1+ I 2wx + — x 

X (l - 2wx - ^x ) 
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we find that the matrix N of (|16|) has the following form 



i 3 / 2 °\ , s 

(pJ+Pc) ! ■ (20) 

\2dp/dp(p c ) v 

Since the eigenvalues of N are independent of the cosmological constant the 
system has a unique bounded solution in the neighbourhood of the centre, which 
is indeed C°° . This implies the existence of a unique, smooth solution to J7J) 
and © or (|T7|) and (|T%|) near the centre. 

In [5] it is claimed that for a fixed equation of state and cosmological constant 
the choice of central pressure (and by virtue of the equation of state the central 
density) do not uniquely determine the solution. The existence and uniqueness 
theorem above disproves this statement. 

Uniqueness of the solution at the centre immediately implies: 

Theorem 2 Let an equation of state p(P), a central pressure P c and the cos- 
mological constant A be given such that 

47rP c +y/ J (P c )-|=0, (21) 

where (47r/3)p c = (47r/3)p(P c ) = w c so that also the central energy density is 
given by the equation of state. Then the unique solution is the Einstein static 
universe with A = Ae- 

For well defined right-hand sides of (|17|) and i|18fl standard theorems for 
differential equations guarantee that the solution can be extended. This im- 
plies that the solution is extendible if the pressure is finite P < oo and if the 
denominator of (|18(1 is non-zero, i.e. y = 1 — 2wx — (A/3)x > 0. Since the sec- 
ond condition depends on the cosmological constant, new properties may arise. 
Moreover it must be clarified whether y = corresponds to a coordinate sin- 
gularity of the spacetime, as in the constant density case 0], or to a geometric 
singularity. In Buchdahl variables H10fl - H12fl the line element JQ) takes the form 

ds 2 = -C(x) 2 dt 2 + — + x(d9 2 + sm 2 6d<j) 2 ) , (22) 

Axy(x) z 

which implies the following non-vanishing components of the Ricmann tensor 

Ret et =R4>t 4>t = -2y(x) 2C -^, (23) 
R x e xe = R x4 > x * = -2y(x)y'(x) , (24) 
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with the remaining two given by 



Re^ = - V{X) \ (25) 



C(i) " w C(*) 
-WW^, (26) 

CO) 

which indicates that y — > corresponds to a coordinate singularity rather then 
a geometric singularity if y'{x) and l/£(x) behave well as y — > 0. 

The following theorem clarifies the extendibility of solutions. 

Theorem 3 Suppose the pressure is decreasing near the centre, which means 

47rP c + ^(P c )-| >0. (27) 

Then the solution is extendible and the pressure is monotonically decreasing if 

AnP + w - A/3 > 0. (28) 

Proof. Assume that p = p(P), P c and A are given such that P is decreasing 
near the centre, then w(x) :X < 0. Since y > equation (fH|) implies 

(tfC,*), x <0. (29) 

Rewriting l|13fl gives 

V(,* = S- (4nP + w-^\ , (30) 



2y V 3, 

next using the implication of (|29|) leads to 

i/C* < (»Cx) («)• (31) 
Together with the explicit expression of y^, x in 130fl this finally shows 

4ttP + w - 4 

2/ > nr- 32 

" 4ttP c + w c - f 

Therefore the Buchdahl variable y cannot vanish before the numerator does and 
consequently the right-hand sides of (|17() and (|18(l are well defined, hence one 
can extend the solution if 4ttP + w — A/3 > 0. 

Since y > and AirP + w — A/3 > the sign of the right-hand side of l|18fl 
is strictly negative. Therefore the energy density and because of the equation 
of state the pressure are monotonically decreasing functions. □ 

The last theorem contained the extendibility condition of solutions, which also 
allows to show the existence of global solutions. 
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Theorem 4 Suppose an equation of state is given such that p is defined for 
p > 0, non-negative and continuous for p > 0, C°° for p > and suppose that 
dp I dp > for p > 0. Furthermore assume that the cosmological constant be 
given such that A < A-Kpi,. 1 

Then the pressure is decreasing near the centre and there exists for any pos- 
itive value of the central pressure P c a unique inextendible static and spherically 
symmetric solution of Einstein's field equations with cosmological constant hav- 
ing a perfect fluid source with equation of state p(P). 

If A < the matter either occupies the whole spacetime with p tending to 
Poo as r tends to infinity or the matter has finite extent. In the second case a 
unique Schwarzschild-anti-de Sitter solution attached as an exterior field. 

If the cosmological constant satisfies < A < kxpb, the matter has always fi- 
nite extent and a unique Schwarzschild-de Sitter solution attached as an exterior 
field. 

Proof. If the cosmological constant is given such that A < Airpb then 

47T A 47T . , A 

< -jPb - - < 4ttP c + y p(iy - 3 , 

and the pressure is decreasing near the centre by J22J). 

Since the pressure is decreasing near the centre the denominator of (|32() is 
some positive number. Furthermore one can estimate the numerator of Ij32|) by 

4ttP + w - 4 
V ~ 47rP c + w c - f ' 

> X - (33) 
" 4ttP c + Wc - f 

47T _ A 

> ,_ n 3 3 a ■ (34) 



4ttP c + w c - § 



and conclude that if 



A < 4irp b , (35) 

the Buchdahl variable y cannot vanish before the pressure does. The coordinate 
Xb where the pressure vanishes will be taken as the definition of the stellar 
object's radius R. 



A < 

If A < the matter can occupy the whole space because (|32(l implies positivity 
of y. 

1 Assumptions on the equation of state could be weakened according to |5] and Moreover 
the line of argument presented in |S] can surely be applied for cosmological constants having 
the derived bound. 
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Suppose that P{xb) = 0, at the corresponding radius R the Schwarzschild- 
anti-de Sitter solution is joined uniquely by the condition M = m(R). In this 
manner the metric is C° only, because the density at the boundary may be 
non-zero. The metric is C 1 at P{R) = if Gauss coordinates relative to the 
hypersurface P(R) = are used. If the boundary density does not vanish the 
Ricci tensor has a discontinuity. Thus the metric is at most C 1 . Without further 
assumptions on the boundary density this cannot be improved. 

Now assume that P(x) > for all x > 0. P(x) is monotonically decreasing, 
therefore ]imx—Kx> P(x) = Poo exists. This implies that the pressure gradient 
tends to zero as x — > oo. Because y~ l — > as x — > oo equation (JJJ does not 
imply that Poo = 0, which it does if A = 0. Thus the equation of state only 
gives p ^ poo = p(Poo) as x — » oo. 

< A < 4np b 

If < A < Airpi, then one can estimate the pressure at the possible coordinate 
or geometric singularity f defined by y(r) = since 



Hence there exists R such that P(R) = 0. Since the pressure is decreasing and 
P(f) < it follows that R<f. 

Thus if the cosmological constant is positive and A < Airpi, then the pressure 
always vanishes at some Xb- At the corresponding radius R the Schwarzschild- 
de Sitter solution is joined uniquely by the same condition M = m(R). The 
metric is at most C 1 because the boundary density is larger than zero, this 
cannot be improved. □ 

It is quite remarkable that this upper bound A < Airpt was found indepen- 
dently by considering consistency of the Newtonian limit [S], the gravitational 
equilibrium via the virial theorem |1U) and was also found in |15) by demand- 
ing stability of circular orbits. This "coincidence" deserves a closer inspection, 
which was already initiated in [Q. 

4 Buchdahl inequality 

The importance of Buchdahl's inequality was already discussed in the introduc- 
tion. In the following its generalisation with cosmological constant is derived, 
see e.g. |3] without the cosmological term and [7] with cosmological constant 
following Buchdahl's original approach. 




(36) 



which holds for all r. Therefore 




(37) 



In the proof of the following theorem the field equations with constant den- 
sity, denoted with a tilde, are solved with the help of Buchdahl variables and 
then this solution is compared with a general decreasing solution. The bound- 
ary mean density of the general solution defines the constant density solution 
that is used for comparison. 

Theorem 5 Let the cosmological constant be given such that A < A-Kpij. Then 
for solutions having finite radius there holds 



Proof. Assume that p(P), P c and A are given such that the pressure is de- 
creasing near the centre. Then by theorem|31the pressure and mean density are 
decreasing functions and equation l|14f) implies 

{K,x),x = , K.x = - ^ , (39) 

where the constant of integration is obtained from H13f) evaluated at the bound- 
ary which can be used to integrate ||2HJ). Let us rewrite the right-hand side of 
dSHJl as 

*•-(»+ 1) G« - 7) (40) 

and substitute in —2yy x = 2w + A/3 for the first factor, divide by y which after 
integration yields 

CW--^-4^)+C(0) + |r-4. (41) 
2w + j 2w + j 

From equation l|29[) we conclude y^ x = {yC,x) b < yC,%> so that with y > y one 
finds 

w-- ~ w- - 

C(*)^C(*) = -^y(x) + C(0) + ^. (42) 

Evaluating this at the boundary and using that £ can be normalised such that 
Cb = Vb, h is found that 

w - A . w- A 

»*-^» + m + 2sTfx- (43) 

Since £(0) is a positive number, some algebra finally leads to the generalised 
Buchdahl inequality 

1 A 

^3-9^ (44) 



Now we compare solutions with decreasing mean density and solution with 
constant density, where the constant density corresponds to the boundary mean 
density of the general solution. Thus (|44|) yields 



1 A 



(45) 



which reads more explicitly 



(46) 



and which holds for all monotonically decreasing densities. 



□ 



By using the definition of the boundary mean density Wb = M/R one can 
reformulate (|46(l arrive at 



a generalised Buchdahl inequality that takes the cosmological constant into 
account. In the 2M/R version of (|47|) was for example used to calculate 
the surface redshift with A. For comparison, note the in [7] the cosmological 
constant is rescaled by 87r. 

5 Solutions without singularities 

The last two sections showed the existence of stellar models for cosmological 
constants A < Airpb- Equation l|47|l implies that the boundary of the stellar 
object has a lower bound given by the black-hole event horizon and an upper 
bound given by the cosmological event horizon. The upper bound for positive 
cosmological constants. 

Stellar models with A < have a lower bound given by the black-hole event 
horizon. At the boundary the Schwarzschild-anti-de Sitter solution attached as 
an exterior field. Therefore stellar models with A < have no singularities. 
Solutions with cosmological constant satisfying A < are globally static. 

For positive cosmological constants the situation is different. But one can 
also construct solutions without singularities. 

At the boundary r = R, defined by P(R) = the Schwarzschild-de Sitter 
solution is joined C 1 by the usual procedure introducing Gauss coordinates. 
The surface r = R can also be found in the vacuum region where the time-like 
Killing vector is past directed. This means that a second stellar object can be 
put in the spacctimc: 

Remarks on finiteness of the radius 

So far it has been shown that given an equation of state, a central pressure and 
a cosmological constant there exist a unique model with finite or infinite extent. 




(47) 
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r = oo, 




r = oo, J? 



Figure 1: Penrose-Carter diagram with two stellar objects having radii R which 
lie between the two horizons. Since the group orbits are increasing up to R the 
vacuum part contains the cosmological event horizon r++. This solution with 
two objects has no singularities. Because of regions lie and II Ic this spacetime 
is not globally static. 

This depends on the given equation of state and on the cosmological constant. 
As was proved in theorem^ solutions are always finite for positive cosmological 
constants. If A < either the pressure vanishes for some finite radius or the 
density is always positive and tends to poo as r tends to infinity, which depends 
on the equation of state; in |12) necessary and sufficient conditions can be found. 

6 Conclusions and Outlook 

In this paper existence and uniqueness of static and spherically symmetric per- 
fect fluid solutions with given equation of state for cosmological constants sat- 
isfying A < Anpb was proven. It would be interesting to extend these proofs 
to any value of the cosmological constant. The existence of these solutions is 
conjectured in £Q. 

One could start studying this system numerically or by using generating 
function techniques. An investigation of the Riemann tensor as outlined might 
also be a good starting point. 

Furthermore the relation of the different theories that give rise to the same 
upper bound must be analysed and understood. 
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